GEOMETRY OF WEAK STABILITY BOUNDARIES 
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Abstract. The notion of a weak stability boundary has been successfully used to 
design low energy trajectories from the Earth to the Moon. The structure of this 
boundary has been investigated in a number of studies, where partial results have 
been obtained. We propose a generalization of the weak stability boundary. We prove 
analytically that, in the context of the planar circular restricted three-body problem, 
under certain conditions on the mass ratio of the primaries and on the energy, the 
weak stability boundary about the heavier primary coincides with a branch of the 
global stable manifold of the Lyapunov orbit about one of the Lagrange points. 



1. Introduction 

We consider the planar circular restricted three-body problem for a small mass ratio 
of the primaries. We give a general definition of the weak stability boundary set in the 
region of the heavier primary. We consider the global stable manifold of the Lyapunov 
orbit about the Lagrange point located between the primaries. We prove analytically 
that, under restrictions on the energy, the weak stability boundary coincides with the 
branch of the global stable manifold in the region of the heavier primary. 

The concept of WSB was introduced in [ll [2] to design low energy transfers from 
Earth to Moon, and subsequently applied to the rescue of the Japanese mission Hiten in 
19910 (See also [3].) A particular feature of the 'WSB method' useful for applications 
is that it allows the capture of a spacecraft into an elliptic orbit about the Moon, with 
specified eccentricity of the ellipse, and with specified true anomaly at the capture. 

There has been considerable work devoted to understand the concept of WSB from 
the point of view of dynamical systems, and to enhance its applicability (see, e.g., [6l[T3l 
m [12]). A remarkable property of the WSB is that, in the context of the planar circular 
restricted three-body problem, for some range of energies, and under some topological 
conditions on the hyperbolic invariant manifolds associated to the libration points, 
the weak stability boundary points coincide with the points on the stable manifolds 
satisfying some additional conditions. This has been observed numerically in [6], and 
argued geometrically in 

The classical definition of the WSB is as follows: for each radial segment emanating 
from the Moon, we consider trajectories that leave that segment at the periapsis of an 
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osculating ellipse whose semi-major axis is a part of the radial segment; a trajectory is 
called weakly n-stable if it makes n full turns around the Moon without going around 
the Earth, and it has negative Kepler energy when it returns to the radial segment; if 
the trajectory is weakly (n — l)-stable but fails to be weakly n-stable, it is called weakly 
n-unstable; the points that make the transition from the weakly n-stable regime to the 
weakly n-unstable regime are by definition the points of the WSB of order n. 

We note that WSB points lie on different Hamiltonian energy levels. Also, the 
WSB is not an invariant set for the Hamiltonian flow. We remark that, since the 
stability/instability criteria, as described above, are concerned with the behavior of 
trajectories for finite time, they inherently introduce 'artifacts', i.e., points with very 
similar trajectories that are categorized differently with respect to these criteria. See 

mm- 

In the present note, we propose a more general definition of the WSB. We remove 
the condition that the infinitesimal mass leaves the radial segment at the periapsis 
of an osculating ellipse whose semi-major axis is a part of the radial segment. We 
remove the condition on negative Kepler energy at the return. We define a point on 
the radial segment as being weakly n-stable provided that it makes n turns around the 
primary, such that the distance from the infinitesimal mass to the primary measured 
along the trajectory does not get bigger than some critical distance. Otherwise the 
point is redeemed as unstable. (Some of these ideas are also suggested in |TT].) The 
main result of this paper is that the WSB points, which make the transition from the 
weakly stable to the weakly unstable regime, are the points on the stable manifold of 
the Lyapunov orbit for the corresponding energy level. 

The argument for the main result is analytical, relying on topological arguments and 
estimates from [5} 110^ [9]. For this reason, we deal with the WSB set about the heavier 
primary (unlike in the WSB original setting). 

An interesting aspect of the WSB method is that it uses 'local' information on the 
dynamics, namely the return of trajectories to a surface of section about one of the 
primaries, to infer some 'global' information on the dynamics, namely the existence of 
trajectories that execute transfers from one primary to the other. 

2. Background 

2.1. The planar circular restricted three-body problem. We consider the planar 
circular restricted three-body problem (PCRTBP) with the mass ratio of the primaries 
sufficiently small. The system consists of two mass points Pi,P2, called primaries, 
of masses mi > m2 > 0, respectively, that move under mutual Newtonian gravity 
on circular orbits about their barycenter, and a third point P3, of infinitesimal mass, 
that moves in the same plane as the primaries under their gravitational influence, but 
without exerting any influence on them. Let // = 1712/ {mi + 7712) be the relative mass 
ratio of m2- In the sequel, we will assume that < ;U < 1 is very small, which will be 
made precise later. 

It is customary to study the motion of the infinitesimal mass in a co-rotating system 
of coordinates (x, y) that rotates with the primaries. Relative to this system, Pi is 
positioned at (//, 0) and P2 is positioned at (—1 + /x, 0). After some rescaling, the 
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equations of motions are given by 

(2.1) x-2y = —, y + 2± = —, 
where the effective potential uj is given by 

(2.2) u;{x, y) = \{x'' + y') + ^ + ^ + ^^(1 " Z^), 

2 ri r2 2 

with n = {{x - iif + y2)i/2^ r2 = {{x + l- fif + y2)V2. 

The equations of motion can be described by a Hamiltonian system given by the 
following Hamiltonian (energy function): 

I 

(2.3) H{x,y,px,Py) = -{{Px+yf + {py-xf) -uj{x,y), 

where x = px + y and y = py — x. 

For each fixed value H of the Hamiltonian, the energy hypersurface Mh is a non- 
compact 3-dimensional manifold in the 4-dimensional phase space. The projection of 
the energy hypersurface onto the configuration space (x, y) is called a Hill's region, and 
its boundary is a zero velocity curve. See Fig. [TJ Every trajectory is confined to the 
Hill's region corresponding to the energy level of that trajectory. 

The equilibrium points of the differential equations (j2.ip are given by the critical 
points of uj. There are five equilibrium points for this problem: three of them, Li, L2 
and L3, are collinear with the primaries (where Li is between L2 and -L3), while the 
other two, L^^L^, form equilateral triangles with the primaries. The distance from Li 
to P2 is given by the only positive solution x+ to Euler's quintic equation 

(2.4) - (3 - ^J)x^ + (3 - 2/i)x^ - ^ix^ + 2fix-fi = 0, 

and so the distance from Li to Pi is 1 — X4-. 

The values H{Li) of the Hamiltonian ()2.3p at the points Lj, i = 1,...,5, satisfy 
iJ(L5) = H{Li) > H{L^) > H{L2) > H{Li). For H < H{Li), the Hill's region has 
three components: two bounded components, one about Pi and the other about P2, 
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and a third component which is unbounded. For H € {H (Li) , H {L2)) , the Hih's region 
has two components, one bounded, which is topologically equivalent to the connected 
sum of the two bounded components from the case H < H(Li), and the other one 
unbounded (Fig. [T]). 

The hnearized stabiUty of the equihbrium point Li is of saddle-center type, with 
the linearized equations possessing a pair of non-zero real eigenvalues itA, and a pair 
of complex conjugate, purely imaginary eigenvalues ziziv. For each H ^ H(Li), near 
the equilibrium point Li there exists a unique hyperbolic periodic orbit 'Jh, referred 
as a Lyapunov orbit. This orbit has 2-dimensional stable and unstable manifolds 
W^{'Jh), W^"(7h), respectively, that are locally diffeomorphic to 2-dimensional cylin- 
ders. These manifolds have the following separatrix property: when restricted to a 
compact neighborhood B}{{cL,b) of 'Jh in the energy hypersurface Mh, of the type 
BH{a,b) = {a < X < b}, with a < x^^ < b sufficiently close to xl-^, each of the 
manifolds W^{'jH),W^{'yH) separates Bn^a^b) into two connected components. 

2.2. Conley's isolating block. Let (/):MxM^'Mbea C-'^-flow on a C^-differentiable 
manifold M. Given a compact submanifold with boundary B C M, with d\m.{B) = 
dim(M), we define 

B- = {pedB\3e>Q s.t. (A(o,£)(p) n B = 0}, 
5+ = {p e 95 I 3e > s.t. 0(_e,o)(p) n B = 0}, 
S° = {p G dB I (pt is tangent to dB at p}. 

We obviously have dB = B^ \J B^ U B^ . We call B~ the exit set and B^ the entry 
set of B. 

An open set V is called an isolating neighborhood for the flow if dV contains no 
orbit of cj). An invariant set S for the flow (p is an isolated invariant set if there exists 
an isolating neighborhood V for the flow such that S is the maximal invariant set in V . 
The compact submanifold B is called an isolating block for the flow (j) provided that: 

(i) B- r\B+ = 5°, 

(ii) B^ is a smooth submanifold of dB of codimension 1, and, consequently, B~ ^ B~^ 
are submanifolds with common boundary B^ . 

The interior of an isolating block is an isolating neighborhood and so determines an 
isolated invariant set, possibly empty. 

In the PCRTBP, Conley has constructed an isolating block around Li that can be 
used to study the nearby dynamics. Consider the part of the Hill's region which satisfies 
a < X < b, where (a, 6) contains the x-coordinate x^^ of Li. This set determines a 
"dynamical channel" which allows for the transit of trajectories between the Pi and 
P2 regions. The lift Bh = BHia,b) of this set to the energy hypersurface, where 
a,b are chosen close to is Conley's isolating block. Geometrically, this is a 3- 
dimensional manifold with boundary dBu consisting of the set of points in the energy 
hypersurface that projects onto x = a and x = 6 in the configuration space. It is 
diffeomorphic to the product of a line segment with a two sphere, Bh ~ [a, b] x 5^, and 
its boundary OBh is diffeomorphic to the union of two 2-spheres, OBh = BH,a^BH,b ~ 
({a} X 52) U ({6} X S^). 
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The isolating block conditions in this case are that every trajectory intersecting dB 
tangentially must lie outside of Bh both before and after the intersection, that is, if 
x{t) = a and x(t) = then x{t) < and if x{t) = b, and x{t) = then x{t) > 0. So we 
have 

Bjj = {{x, y, X, ij) S BBh \ x{t) = a or x{t) = h and x{t) = 0}, 

BJ^ = {(x, y, X, y) E OBh j x{t) = a and x{t) < 0, or x{t) = b and x(t) > 0}, 

B^ = {{x, y, X, y) G BBh \ x{t) = a and x{t) > 0, or x{t) = b and x(t) < 0}. 

For each component of dBj^, the exit and entry sets determine a pair of disjoint 
open 2-dimensional topological disks, which we denote as follows: B^^, BJj^ are the 
exit sets of the boundary components Bn^a, Bjj i,, respectively, and B^ ^, B^^ are the 
entry sets of the boundary components Bn^a, Bh^^ respectively. The complement in 
Bhi, of -B^^ U B^^ is the set -B^^ = -B^ fl {x = 6}. A similar statement holds for 

BH,a- 

The exit and entry sets are further broken up into components with dynamical 
roles. The set B~^^ is the union of three sets, a spherical cap -B]^'^, corresponding to 
trajectories that enter the block Bh through the entry part of -B_f/,b and later leave the 
block through the exit part of B^ a, a spherical zone -B^'^, corresponding to trajectories 
that enter the block Bh through the entry part of Bn^b aiid leave the block through 
the exit part of Bn^b^ and a topological circle separating them, corresponding to the 
intersection of VF^(7i?) with Bufi- Similarly, BJj ^ = B^'^ U B^^ U {BH,b n VF"(7i?)), 
where the notation is analogous to the above. There is a similar decomposition for the 
entry and exit set components of BH,a- See Fig. [2l 

Later in the paper, we will use the following fact, which is a consequence of the 
above discussion. There are three possible behaviors for trajectories that start from 
the Pi-region and enter the isolating block: 

(i) Trajectories enter the block through -B^'^, exit the block through B^'^^, and so 
they execute a transfer from the Pi-region to the P2-region. 

(ii) Trajectories enter the block through -B^'^, exit the block through -B^'^, and so 
they do not transfer to the P2-region. 

(iii) Trajectories enter the block through -B//,6niy (7//) and are forward asymptotic 
to 7h, and so they never leave the block. 

For further details on this subsection, see [5]- 

2.3. Hyperbolic invariant manifolds. The geometry of the hyperbolic invariant 
manifolds can be described analytically inside the Pi-region, for some range of energies 
and mass ratios, following some results from jlO^ [9]. 

First, there exists an open set Oi in the (/i, i7)-parameter plane, with < /x <C 1 
and H ^ H{Li) such that, for {fJ-,H) € Oi, the following hold: 

(i) The energy hypersurface Mh contains an invariant 2-torus Th separating Pi 
from Li. 

(ii) There exist a < xli < b such that the flow inside the isolating block Bh = 
BH{a,b) is conjugate to the linearized flow. 
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(iii) In the region Nh in Mh bounded by Th and Bn^b-, the longitudinal angular 
coordinate 6 is increasing along trajectories. 

Second, for all < /x ^ 1 sufficiently small, the (x, y)-projections of the branches of 
Ty^(-Li) inside the Pi-region have the following properties: 

(iv) The distance d to the zero velocity curve, and the angular coordinate 9 satisfy 
the following estimates: 

(2.5) d = n^'^ (fiV - 3^/6 + M cos t + o(l)) , 

(2.6) 9 = --K + n'^/'^ {Nt + 2Ms\nt + o{l)) , 

where M, are constants, the parameter t means the physical time measured 
from a suitable origin, and o(l) — > when ^ — )■ uniformly in t as i = 0(/i~^/^). 
These expressions hold true outside Bh- 

(v) There exists an open set O2 ^ Oi in the {jj,, i?)-parameter plane, with < 

<C 1 and H ^ H{Li) such that, for {n,H) G O2, the (x, y)-projections of 
the branches of W^^jh), W^^jh) inside the Pi-region satisfy estimates similar 
to (j2.5p and (j2.6p . That is, these invariant manifolds turn around Pi in the 
region Mh bounded by the torus Th and the boundary component -B//,6 of the 
isolating block Bh- Moreover, there exists a sequence of mass ratios for 
which W'^{'~fH) and W^{'~fH) have symmetric transverse intersections, provided 
{fik,H) G 02- 

The geometry of the hyperbolic invariant manifolds for the range of parameters 
considered above allows to extend the separatrix property of these manifolds from the 
local clS described in Subsection 12.11 to the global case. For as long as the stable 
and unstable manifolds do not intersect each other, the cuts of these manifolds with 
a surface of section are topological circles. If a point is inside the i-th cut Tg^ {(ih) 
made by the stable manifold W^^(7_f/) with the surface of section Sq^, which is assumed 
to be a topological circle, then the forward trajectory of that point stays inside the 
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Figure 3. Projection of McGehee's separating torus onto configuration 
space, and trajectory near the zero velocity curve. 



cylinder bounded by W^{^h) in Mh for z-turns and transfers from the Pi-region to the 
P2-i'egion afterwards. If a point in Sg^ is outside the i-th. cut Tg^ iiln)-: then its forward 
trajectory stays inside the Pi-region for at least (z + l)-turns. A similar statement holds 
for the cuts made by the unstable manifold and backwards trajectories. 

If the stable and unstable manifolds intersect, say j(7H) intersects r^(jj(7//)j 
then the intersection points are homoclinic points that make {i + j)-turns about Pi, 
and some future cuts of the invariant manifolds cease to be topological circles. For 
example, Tq^^j^-{'^h) is a finite union of open curve segments whose endpoints wind 
asymptotically towards F^^ ^(7//). Due to the asymptotic behavior of the endpoints, 
each of these open curves divides Sq^^ into transfer and non-transfer orbits. Thus, 
the separatrix property extends to the case when the cuts of the hyperbolic invariant 
manifolds cease being topological circles. See OH]. 

There are no analogues of the above analytical results for the P2-region about the 
lighter mass. 



2.4. Equations of motion relative to polar coordinates. We recall the relations 
between the motion of the infinitesimal mass P3 relative relative to the barycentric 
rotating coordinates {x,y,x,y), relative to the polar coordinates {r,9,r,6) about Pi, 
and relative to the classical orbital elements {a,e,(l),T) about Pi. 

The relation between barycentric and polar coordinate is r = ((x — /i)^ + y^)^^^ and 
tan^ = y/{x — 

The orbital elements are characterized by the semi- major axis a of an ellipse with a 
focus at Pi, the ellipse eccentricity e € [0, 1), the argument of the periapsisis cp £ [0, 2tt], 
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and the true anomaly r E [0, 27r]. We have the following coordinate transformations 
X = r cos((/) + r) + //, 
y = r sm(6 + r), 

(2.7) 

i; = r cos(0 + r) — rf sin((/) + r) + r sin(i?i> + r) , 
y = r sin((/) + r) + rf cos((/) + r) — r cos(i?^ + r), 
and the following formulas 

a(l - e^) 

r = 

1 + e cos r ' 

ae(l — e^)fsinT 

(2.8) ^ ~ (l + ecosT)2 ' 

= + T, 

. . Va(l-e2)(l-/i) 

U = T = 2 • 



The Hamiltonian function in polar coordinates is given by 
2.9) H{r,pr,e,pg) = ;^(p^ + \e^) -pe + f^r cos 6 + uj{r,pr, 9, pe) 

where the canonical momenta are given by 

Pr = r,pe = r^iO + 1). 

Note that the conservation of energy implies that the initial position (r, 9) relative 
to Pi and the initial radial velocity r uniquely determine a trajectory, up to a choice 
of a sign for 6. Suppose that we know the initial data (r, r, 9) on a trajectory. Using 
(|2.8p . the eccentricity of the osculating ellipse to this trajectory at the initial point 
uniquely determines the trajectory, and hence its energy. This implicitly defines (p and 
T. Conversely, if we have a trajectory for which the initial angle coordinate 9, the 
initial angular velocity r, and the eccentricity of the osculating ellipse e at the initial 
condition are fixed, then the energy level H of the trajectory uniquely determines its 
initial value of r. 

3. Weak Stability Boundary 

We consider the system of polar coordinates (r, 9) about Pi as above, and we let 
H(r, r, 9, 9) be the Hamiltonian relative to this coordinate system. As discussed above, 
the energy is also uniquely determined by the (r, r, 9, e)-data, where e is the eccentricity 
of the osculating ellipse at the initial point. We consider a Poincare section through 
Pi that makes an angle with the x-axis, which is given by 

Se, = {{r,f,9,9)\9 = 9Q,9>{)}. 

Let denote the radial segment obtained as the intersection of Sq^^ with the (x, y)- 
space. Any trajectory that meets Sq^^ transversally is uniquely determined by the (r, r)- 
coordinates of the intersection point, as the ^-coordinate equals in this section, and 
the ^-coordinate can be solved uniquely from the energy condition if(r, r, 9, 9) = H, 
provided 9 > 0. 
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Consider a trajectory with the initial condition zq = ZQ{rQ,rQ,6Q,eo) with initial 
position r(0) = rg, ^(0) = 0q, initial radial velocity r(0) = f-Q, and 9(0) > 0, for 
which the osculating ellipse at the initial point has eccentricity cq- We keep the values 
of rQ,9o,eQ fixed and investigate the change of behavior of the trajectories when ro 
changes. Note that different initial values of rg yield different energies Hq. 

Fix a value of n sufficiently small for which there exists an open range of energies 
{H{Li),H*) with {n,H) £ O2 for each H € {H{Li),H*), as in Subsection ESI For 
this range of energies the estimates (|2.5p are valid. Fix a < x^^ < b such that Bh{cl, b) 
is an isolating block for all H € {H{Li),H*). Let y^, be the supremum of the y- 
coordinates on the segment x = b inside the Hill's regions for H G {H{Li),H*). Define 
61 = arctan(y;,/(^ ~ b))- Let Di be the distance from Pi to x = a, that is Di = ^ — a. 

Fix H € {H{Li), H*) and consider the projection V'^(x,y){-^H) of Mh onto the (x, y)- 
configuration plane. For each angle coordinate 6 G [0,27r], there exists a well de- 
fined interval {ri{H,6),r2{H,9)) such that {r,9) G P^{x,y){-^H) if and only if r G 
{ri{H,9),r2{H,9)). For each trajectory point {r,9) G P^(x,y){-^H) there exists a set 
of admissible values of the radial velocity r and of the eccentricity of the osculat- 
ing ellipse e corresponding to the trajectory at that point. When we let H vary in 
{H{Li), H*), then for each 9 G [0,27r], we obtain an open set of admissible values of 
(r, e) corresponding to all trajectories for all of these energy levels. 

We fix an angle ^0 and a pair of admissible values {ro,eo). Since the energy H is 
uniquely determined by the data (ro, vq, 9o, cq), there exists an open set TZ{ro, 9o,eo) C 
{ri{H,9),r2{H,9)) of ro-values such that //q = H{rQ,rQ,9Q,eo) G {H{Li), H*) provided 
ro G TZ{'rQ,9o, eo). In the next definition, we will consider trajectories with initial points 
Zq lying on the radial segment Iq^ . We will restrict to values of ro in the set 7^(ro > ^0 > ^o ) • 

Definition 3.1. We say that a forward trajectory with initial point zq = zo{ro,9Q) in 
, initial radial velocity ro and initial eccentricity of the osculating ellipse eo , is weakly 
n-stable provided that it turns n-times around Pi , with all intersections with Ig^ being 
transverse, and such that the distance to Pi is always less than Di. If the trajectory 
is weakly (n — l)-stable but fails to be weakly n-stable, we say that the trajectory is 
weakly n-unstable. 

The conditions on the parameters assumed for the Definition 13.11 are imposed in 
order to define the critical distance Di in a consistent way for the whole range of 
energy values H G {H{Li), H*). We recall that in the classical definition of the WSB, 
a trajectory is called n-stable if it turns n-times around Pi, without turning around 
P2; in that case one can consider the distance from Pi to P2 as the critical distance. 

We note that the transversality requirement in Definition 13.11 on the intersections 
of the trajectory of the infinitesimal mass with Iq^, implies that weak n-stability is an 
open condition, that is, if a trajectory starting at some zq = {ro,fo,9o,eo) is weakly 
n-stable, then all trajectory starting inside some domain of the type 

(r, r, 9, e) G (ro - e, ro + e) x (ro - e, ro + e) x (6*0 -e,9o + e) x (eo -e,eo + e) 

with e > sufficiently small, are also weakly n-stable. 

Thus we obtain the following set of weakly stable points in the phase space 

Wn = {zo{ro,rQ,9o,eo) \ zq is weakly n-stable relative to l0^^, 9o £ [0, 2tt]}. 
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Due to the open conditions on the n-stable trajectories, the set Wn is an open set 
of points in the phase space. If we fix the parameters ro, Oq and eo, then we obtain an 
open set Wn(ro, 6*0) cq) in Zg,,, which is a countable union of disjoint open intervals 



The points of the type r2k-i,'f2k at the ends of these intervals are weakly n-unstable. 

Definition 3.2. The WSB of order n, denoted W*, is the set of all points r*(ro, ro; ^O; ^o) 
that are at the boundary of the set of the weakly n-stable points, i.e.. 



We also denote by W*(roi ^o^ ^o) the set of WSB points on the radial segment Ig^ of 
fixed parameters ro and eo- Thus, the WSB set W*(ro;^0;eo) contains the closure of 
the set of all points of the type r2k-i,i"2k, which are the endpoints of the intervals of 
weakly n-stable points within each radial segment as in (13. ip . 

The main result of the paper says that, if we restrict to some angle range of 
outside the angle sector [n — 9i,tt + 9i], where 6i is defined as above, then the WSB 
set is completely determined by the stable manifolds of Lyapunov orbits. To state 
this result, we have to adopt a convention on how to count the number of cuts made 
by the stable manifold with a surface of section Sg^. We label a cut made by the 
stable manifold W'^{'yHo) with Sg^ as the i-th. cut provided that the net change AO of 
the angle 9 along all trajectories starting from Sg^ and ending asymptotically at 'JEq 
satisfies 2i7r < A9 < 2{i + l)ir. Note that as long as Oq ^ [tt — 0i,7r + ^i] there is no 
ambiguity about the labeling of the cuts with the section Sg^^ . 

Theorem 3.3. Fix a pair of admissible values (ro,eo) as defined above. Assume 6q G 
{—IT + 0i,TT — 6i), where 9i is defined as above. Then a point zq = ZQ{rQ,ro,6Q,eo), 
with ro € TZ{rQ,9o,eo), is in yV*{ro,0Q,eo) if and only if zq lies on the (n — l)-st cut 
n-iilHo) of the stable manifold W'^{'yHo) with the surface of section Sg^, where Hq 
is the energy level corresponding to zq. 

The restrictions imposed on the parameters in Theorem 13.31 are needed to apply 
the analytical arguments from Subsection 12.31 It is nevertheless shown in [3] that the 
WSB overlaps with some subset of the stable manifold of the Lyapunov orbit under 
much weaker conditions, provided that the hyperbolic invariant manifolds satisfy some 
topological condition (they turn around the primaries for a long enough time, without 
colliding with the primaries). Moreover, in [3] a wider energy range is considered, in 
which case the WSB is identified with a subset of the union of the stable manifolds of 
the Lyapunov orbits about Li and about L2. The situation described by Theorem l3.3l is 
just a special case when the required topological conditions can be verified analytically. 

Now we explain the relation between WSB and hyperbolic invariant manifolds in 
a more concrete way. Assume that we fix some energy level Hq G {H{Li),H*). We 
generate the stable manifold of the Lyapunov orbit ^Hq: and we count the successive 
cuts made by the stable manifold with some Poincare surface of section Sg^^ . Let zq be 
a point on the (n — l)-st cut Tg^ „_i(7_h'o) of W'^{'yHo) with Sg^^. Let ro be the radial 
velocity at zq, and eo the eccentricity of the osculating ellipse at zq. Then the point zq 



(3.1) 




k>l 



W: = dWn. 
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is in the WSB set >V*(ro> ^O) eo)- Moreover, every WSB point can be obtained in this 
way. 

4. Proof of the main result 

Due to the angle restriction 6*0 G (— vr + 6'i,7r — ^i), in Theorem \3.3\ we restrict to 
the following set of weakly n-stable points 

y^n = {-20 ('^0) ''^0) ^Oi Co) I is Weakly n-stable relative to Ig^, 9q € (— vr + ^i, tt — 6*1)}. 

Since the n-stability is an open condition and the angle range (— vr + 0i,7r — 9i) is 
also open, the set Wn is an open set in the phase space. 

We prove that a point zq = zo{ro,ro, Bq, cq) is in W* if and only if it is in the (n — 1)- 
st cut Tg^ 1(7^0) made by the stable manifold W^^jho) with Sqq, where Hq is the 
energy level corresponding to zo- For this, we first show that zq is a weakly n-stable 
point on Ig^^ if and only if it is outside the domain in Sg^ bounded by Fg^ „_i(7j?o)' is 
weakly n-unstable if and only if it is inside the domain in Sg^^ bounded by Fg^ „_i(7_ffo)- 

First, we show that the points inside the cylinder bounded by the stable manifold 
are weakly unstable. Let zq = zo{ro,rQ,9o,eo) be a point in /gg. Then (j2.9p gives 
the value Hq of the energy of the trajectory with initial condition zq. Assume that 
Zq is inside the domain in Sg^^ bounded by Fg^ „_i(7_h'o)- the separatrix property 
from Subsection 12.31 the trajectory turns counterclockwise precisely (n — l)-times inside 
the domain Mho, while staying inside the region of the cylinder bounded by W'^{^Ho)j 
enters the isolating block through the entry set region B^'^^^, crosses the block 

and exits it through the exit set region BJ^'^^. When the trajectory leaves the block 
Bhq, the distance from Pi is bigger than Di. Since the trajectory achieves a distance 
to Pi bigger than the threshold value prior to completing an n-th turn around Pi, the 
trajectory is weakly n-unstable. 

Second, we show that the points outside the cylinder bounded by the stable man- 
ifold are weakly stable. Assume that zq is outside the domain in Sg^ bounded by 
Fgg „_i(7Ho)- the separatrix property from Subsection 12.31 the trajectory will turn 
counterclockwise inside the domain Mhq and will keep staying outside the region of the 
cylinder bounded by W'^{^Ho) for at least n turns. If the trajectory leaves the domain 
Mhq, it has to meet the block Bh^ at or at P^^'^ ^. In the first case, the trajectory 
bounces back to the domain Mhq and it continues its counterclockwise motion about 
Pi. In the second case, it cannot leave the block Bh^ through BHo,a, since only the 
points that are inside the cylinder bounded by W^{'yHo) can do that; it cannot remain 
inside the block Bhq for all future times since only the points on W^^jho) have this 
property; hence, it has to leave Bh^ through the exit set region Pj^^^f,, and to go back to 
the domain Mhq- The time spent by the trajectory inside the block Bhq does not affect 
the count of turns about Pi. Since zq is outside the cut Fg^ „_ 1(7^0)' the trajectory 
cannot leave the Pi-region after only (n — l)-turns, so it turns around Pi for at least 
n-turns. Thus the trajectory is weakly n-stable. 

Now, we prove the statement of the main theorem. 



12 



E. BELBRUNOt, M. GIDEA*, AND F. TOPPUTO 



First, assume that zq is in Tg^ „_i(7Ho)- Its forward trajectory turns (n — l)-times 
around Pi and then approaches asymptotically 7^0 • Thus the trajectory is weakly 
n-unstable. To show that zq is an WSB point it is sufficient to prove that there exists 
a sequence {zQ)k>i with Zq G Wn and Zq ^ zq as k ^ 00. Take a small 4-dimensional 
open ball U around zq in the phase space. Let T > be such that the time-T map (pT of 
the Hamiltonian flow takes zq to a point in B}{o,b, where Hq = H{zq). The image (pxipf) 
oiU by (t)T is a 4-dimensional open topological ball about (t)T{zo). We intersect (pri^^) 
with the 4-dimensional submanifold with boundary U_f/G(H(Li) H*) ^H- The ball (prii^) 
has non-empty intersection with [Jh£(h{Li) h*) ^ti\- These intersection points yield 
weakly n-stable trajectories. Thus (pT{U)^\jHe(H{Li) H*) contains a 4-dimensional 
open, topological ball V, which contains 4)t[zq) on its boundary, consisting of points 
that correspond to weakly n-stable trajectories, i.e., those trajectories that return to 
the Pi region for at least one extra turn about Pi. 

Now consider the set 0_r(V). This is a 4-dimensional open, topological ball in lA 
that contains zq on its boundary. There exist 9' arbitrarily close to 6q such that the 
intersection (/)_t(V) fl Sqi is a non-empty open set. All points z' G (p-xiV) fl Sqi are 
weakly n-stable points. We note that these points may not lie on Iq^, nor on the same 
energy level as zq ; they can also have the eccentricity of the osculating ellipse different 
from cq. Thus, arbitrarily near zq one can always find weakly n-stable points, and since 
Zq itself is weakly n-unstable, it follows that zq S dWn = VV*. 

Second, assume that zq G >V*(rO) ^O; cq). Then there exists a sequence of points 
(•2o)fc>i on 1{9q) such that Zq is weakly n-stable and — )• as — )• 00. From the 
above, we know that the weakly n-stable points are those inside the cylinder bounded 
by the stable manifold. Thus, there exists a corresponding sequence of stable manifold 
cuts Fg^ „_i(7_H'^.) where i/^ = H{zq), such that Zq is inside the region in Se^ bounded 
n~i('yHk)- Since z^ — )■ zq it follows that H{zq) H{zo) = Hq as k ^ 00. 
The stable manifold cuts also depend continuously on the energy, so ^0^^ n-ii'^Hk) 
approaches F^^_^_-^(7j^J as /c 00. Hence zq G F^^ ,^_-^(7Ho)• 
Through double inclusion, we conclude that 

Zq G >V*(ro, Oo,eo) if and only if zq G F^^ 1(7^0 )• 

5. Concluding remarks 

We compare the invariant manifold method with the WSB method. The invariant 
manifold method is based on identifying geometric objects that serve as building blocks 
that organize the global dynamics: equilibrium points, periodic orbits, and their stable 
and unstable invariant manifolds, if they exist. The WSB method is a local method 
for deciding whether the trajectories about one of the primaries exhibit some kind of 
stability in terms of the return to a surface of section. The conclusion of this paper, 
corroborated with the results in [6l|4], is that in simple models the two methods overlap 
for a substantial range of parameters. 

One can think of some other kinds of indicators that mark the passage between the 
weakly n-stable and the weakly n-unstable regimes. One such a possible indicator is 
the continuity of the Poincare return map. The n-th return map to Sqq is continuous 
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at all weakly n-stable points. At the WSB points, the return map exhibits essential 
discontinuities of infinite type. Thus, the set of points where the return map fails to 
be continuous contains the WSB points. 
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